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' , I will briefly review the progress which has been made in the investigation of secondary 

heavy quark production at LEP. I will show why a calculation of the secondary heavy 
(— I I quark multiplicity, keeping the dependence on an event shape variable, is essential for better 

exploitation of data, and I will present the results of such a calculation. This will be compared 
to Monte Carlo studies. 

X 

' 1 Introduction 

Heavy quark production in e''~e~ annihilation can come from two possible sources: from the 
hard interaction itself, e+e~ — > QQ, and from the splitting of perturbatively produced gluons, 
e''"e~ — > qqg — > qqQQ. I shall refer to the former as primary heavy quarks and the latter as 
secondary heavy quarks. 

Of course, in order to be able to regard primary and secondary heavy quark production as 
separate processes, the interference between them must be zero, or at least very small. Fortu- 
nately this is the case[Q]. For non- identical quarks coupled to a vector current, the interference 
between them will vanish by Furry's theorem if we assume that the charges of the quarks are 
not measured. For an axial current, cancellations will occur between up and down type quarks, 
leaving only the case where the "light" quark is a bottom quark. This will only provide an 
effect of the order of 0.2% of the secondary heavy quark rateQ. The contribution for identical 
quarks is only slightly larger. In this case Furry's theorem no longer applies but necessarily 
the quarks will be of the same flavour (bottom or charm) and subsequently give only a small 
contribution. Therefore one can regard secondary heavy quark production as a separate process 
and investigate it both experimentally and theoretically. 

In this talk, I will briefly review the progress which has been made in this investigation. I 
will show why a calculation of the secondary heavy quark multiplicity, keeping the dependence 



on an event shape variable, is essential for better exploitation of data, and I will present the 
results of such a calculation. This will be compared to Monte Carlo studies. 

1.1 Theoretical Progress 

The heavy quark mass provides a natural infrared cut-off for the theoretical calculation of the 
total rate of secondary heavy quark production. It is therefore an infrared-safe quantity, and 
can be calculated as an order- by-order perturbative expansion in qs, starting at 0{ag). At 
higher orders in as, large logarithms of s/niq arise, potentially spoiling the convergence of the 
perturbative series at high energies, s ^ m-g. It is interesting to note that these are exactly 
the same logarithms which occur in the resummation of jet multiplicities, where they appear as 
logarithms of the jet resolution scale, ycut- Here we are able to study their effect without the 
arbitrariness of jet algorithms. 

In Ref. the leading and next-to-leading logarithms were summed to all orders in as 
yielding a result that is uniformly reliable for all s. This gave the fraction of decays that 
contain a secondary charm or bottom quark pair, of 

/e = 1.349%, /b = 0.177%. (1) 

1.2 Experimental Progress 

Several experimental measurements of the secondary heavy quark production rate have now 
been made. In Ref. Q, it was extracted for charm quarks from a measurement of the D* 
fragmentation function, and found to be more than a factor of two above the expectation of 
Ref. Q, although with large systematic errors coming from uncertainty in the fragmentation 
function of primary charm quarks. Refs. made less model dependent measurements by 

selecting hard three-jet events, which enhances the fraction of heavy quarks produced by the 
gluon splitting mechanism. In general the measurements have been above the predictions of 
Eq. (|l|), although within the range allowed by variations in as and the quark mass. 
The combined LEP value^ obtained are, 

/c = (2.44 ± 0.43)%, /b = (0.22 ± 0.13)%. (2) 

Of course, in order to extract these results one must use some theoretical input in order 
to separate the primary and secondary heavy quark contributions. In the above experiments 
this theoretical input was the shape of the secondary heavy quark multiplicity distribution with 
respect to some event shape variable. For example, Ref. |^] uses the shape of the multiplicity 
distribution with respect to the jet mass difference. This is fitted to data, extracting a value for 
the overall normalisation and therefore the total rate. Since a full perturbative QCD calculation 
of this shape was not available, Monte Carlo event generators ||7|-^ were used for this theoretical 
input. 

1.3 A New Calculation 

Clearly, it is desirable to calculate the resummed multiplicity, retaining the dependence on the 
jet kinematics, and thereby allowing the calculation of any event shape variable to be performed 
numerically. This perturbative QCD calculation can be made accurate up to next-to-leading 
logarithms whereas the Monte Carlo event generators are only accurate to leading logarithms. 
Furthermore, the leading order, ©(a^), can be included exactly, and the exact kinematics can 
be used when calculating the required event shape. 

"For the bottom quark there has been only one measurement , while for the charm quark, we have averaged 
the results of Refs. flB and 0], assuming that the systematic errors are uncorrelated. 



Here I present the multiphcity of secondary heavy quarks as a function of the heavy jet mass. 
I employ the 'thrust-like' definition of heavy jet mass, where the thrust axis is used to divide 
each event into two hemispheres and the mass of the heavier hemisphere is defined as the heavy 
jet mass. Clearly, the calculation can be easily modified to any thrust-like event variable. 



2 The Secondary Heavy Quark Multiplicity as a Function of Heavy Jet Mass 

The leading order differential cross section for the production of secondary heavy quarks, 
7* — > qqQQ, is easily calculated, and will not be discussed further here. For the calculation 
of the logarithmic contribution, there are two requirements to which we must conform. Firstly 
we must retain the exact kinematics of the qqg production in order to be able to accurately 
obtain the heavy jet mass for each event. We need not worry about the exact kinematics of the 
heavy quarks, since the large logarithms arise from the parts of phase space where they become 
collinear and only the kinematics of the gluon need be considered (although the exact kinemat- 
ics of the heavy quarks are, of course, included in the fixed order contribution). Secondly, we 
must include all soft gluon emission from the light quarks and virtual gluon. These emissions 
contribute to the heavy jet mass (by taking the light quarks off mass-shell) and also provide 
large logarithms which must be summed to all orders using the coherent branching formalism. 
Bearing these considerations in mind, the differential multiplicity is taken to be, 

2 _j_ 2 

n^_(M|f, Q'^;Qq) = J dxi dX2 dkf dk^ dkg ^2 yQ2^A(^l' ^lmax)fqi^h ^2 max) 

n'^^ikl kl;Ql)e{kl - Ql)5{Ml - /i(xi, X2, kl kl fc^)). (3) 

As usual, xi and X2 are the energy fractions of the light quark and antiquark respectively, and ki, 
k2 and kg are the four- momenta of the quarks and gluon. The maximum value of kf , i = 1,2, as 
constrained by the phase space limits, is kf^^^. Also, k']_ is the transverse momentum (squared) 
of the virtual gluon, given by, 

kl = {l-xi+ei- e2)(l - X2 + €2 - ei)Q^, (4) 

where is the centre-of-mass energy squared and e^, (i = 1,2) are the rescaled (primary) quark 
masses (squared) which result from the soft gluon emission. By retaining the exact kinematics of 
the quark, antiquark and gluon, we are able to calculate the heavy jet mass /i(x]^ , X2 , , A^2 5 ^0 ) 
exactly, thereby satisfying the first of our requirements. 

The second requirement, that we should include all soft gluon emission, has been achieved 
by the inclusion of the functions fq and n^^. The function fq is the quark jet mass distribution 
which has been calculated to next-to-leading logarithmic accuracy in Ref. |1C|. More explicitly. 



/g(/c^, Q'^) dk'^ is the probability that a quark, created at a scale gives rise to a jet with mass 
squared between k'^ and k'^ + dk'^ . This function includes all soft gluon emission from the light 
quarks and sums, to all orders, leading and next-to- leading logarithms of kf^^^^/kf, i = 1,2. 
The function n'^^ is the gluon jet mass distribution weighted by the heavy quark pair 

multiplicity. In other words, n'^'^ {kg , k'^; Qq) dkg is the number of heavy quark pairs within 
a gluon jet which was formed at a scale k'j_ and has a mass between kg and kg + dkg. To the 
accuracy required here this quantity must include leading and next-to-leading logarithms of both 
kl/k^ and kl/Ql 

One important simplification of the calculation is allowed by the introduction of a heavy 
quark effective massQ, 

m*Q = imQe^/s. (5) 



Since the bottom quark mass cannot be neglected we should include all mass effects in the 
g ^ bb splitting. This would lead to integrals of the from, 




where x is the rescaled heavy quark mass, mq /-y/s, and x* is the rescaled effective mass. However, 
exactly the same result is obtained (to next-to-leading logarithmic accuracy) if the massless 
splitting function and the effective mass are used. 

lim / — log"^-^ z = -- log" X*' 

Therefore by using this effective mass one can neglect the heavy quark mass in the decay of the 
gluon while maintaining the correct leading and next-to-leading logarithms. In Eq. (P) this is 
manifest as the resolution scale at which the heavy quarks are resolved, Qq = 2mQ. 



(7) 



2.1 The Multiplicity Weighted Mass Distribution 

It is more convenient to calculate the integrated distribution 

10 



N^^{k\ Ql) = £^ dq'nfiq^ Q\ QI). (8) 



Physically this is the number of QQ pairs resolved in gluon jets of mass squared less than fc^. It 
can be derived from N^{k'^, Q^; Qq), the multiplicity of gluons within gluon jets of mass squared 
less than fc^, by integrating over the kernel for the splitting g QQ. Since we have introduced 
the effective mass mg, the appropriate kernel is the massless splitting kernel Pqg. has been 
derived in Ref. |11|. The integration yields, 

NQ<^{e, g2) ^ ^^(^2^ I _^Q<9(^2. + Ca {I'^^ik'; Ql) - l^^ikVQ'; Ql)) } , (9) 

with, 

3bCA I V 4 ^oZk 

-^(2S-2-C)-l(C + 2)log f|)-jfl-| 

In the above, M^'^{k'^]Ql) is the multiplicity of heavy quarks pairs found in a gluon jet, re- 
gardless of the jet mass, and Fg{k'^,Q'^) is the probability that a gluon jet produced at a scale 
Q will have a mass squared below k"^. The notation used above is standard and can be found in 




Ref. |12|. A/i is defined as N but with B replaced by — 1. 

Notice that Eq. (|9|) conforms with nai've expectations. One might expect that the number 
of heavy quark pairs found within a gluon jet created at a scale Q^ and of mass below k"^ would 
be simply the probability of finding a gluon jet of mass below k"^ multiplied by the number of 
heavy quark pairs within. Indeed, this is the first term of our expression for N^^ , and the nai've 
expectation requires only the addition of next-to-leading logarithmic corrections. 

The expression for can now be trivially obtained from the differentiation of A'^^'^ with 
respect to the jet mass. This result is then matched to the fixed order result to avoid double 
counting. 



2.2 Calculation of the Background 



The background to secondary heavy quark production in e"'"e~ annihilation, i.e. primary heavy 
quark production, is estimated by standard three-jet production, since the mass effects will be 
small. Care has been taken to include both the fixed order contribution and all large leading and 
next-to-leading logarithms. As for the secondary heavy quarks, the logarithmic contribution is 
matched to fixed order to avoid double counting. The full result can be seen plotted in Fig. ||. 



3 Numerical Results 



For all the distributions we show, we concentrate on their shape, normalised to the number of 
secondary heavy quarks, rather than on the total rate. We use as = 0.118 and values of 1.2 and 
5 GeV for the charm and bottom quark masses respectively. 

We present the heavy jet mass distribution for \/s = mz in Fig. ||. This is closely related to 
the jet mass difference, Mh — Ml, which was the event shape used to fit fc in Ref. We see 
that the heavy jet mass provides a good discriminator of events with secondary heavy quarks 
from the three-jet background. 
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Figure 1: The multiplicity of primary and secondary Figure 2: The multiplicity of secondary bottom quark 

heavy quark pairs as a function of the heavy jet mass, pairs as a function of the heavy jet mass, normalised 

normalised to the number of heavy quarks. to the number of TP decays. 



Of course, these shapes are dependent on the values chosen for the parameters, Aqcd and 
mq. The effect of varying the quark mass by 5% (dotted) and Aqcd by a factor of two (i.e. 
as by 10%) (dash-dotted) is compared to the result with as = 0.118 and mb = 5 GeV (solid). 
Also shown is the contribution from the fixed order term alone (dashed), demonstrating the 
importance of resumming large logarithms. 



4 Event Generators 



The predictions from Monte Carlo event generators for ^Js = mz can be seen in Fig. ^. Since 
these models are only formally accurate to leading logarithms and do not include the exact 
matrix elements for qqQQ production, our calculation is more accurate and can be used to 
check them. We see that HERWIG and JETSET give similar predictions for the distribution as 
well as for the rate and that ARIADNE peaks at a somewhat lower heavy jet mass. Following 
a suggestion of Ref. later versions of ARIADNE contain an option to veto gluon splitting 
with rUg > k±g. Adding this modification, ARIADNE's distribution is more like the other 
models', but still somewhat different, particularly at low jet masses. Our results lie between the 
unmodified ARIADNE and the other models. 
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Figure 3: The multiplicity of bottom quark pairs as 
a function of the heavy jet mass at -Js = mz- 



Figure 4: The multiplicity of bottom quark pairs as 
a function of the heavy jet mass at = 500 GeV. 



Increasing the centre-of-mass energy, the relative importance of the fixed order term is 
reduced and one gets a cleaner probe of the parton evolution. In Fig. |^ we show the comparison 
at -^s = 500 GeV. The modified version of ARIADNE is in even better agreement with the 
other two models, while the unmodified version is in good agreement with our calculation. We 
therefore see no evidence to support the claim of Ref. [j^ that there is a problem with ARIADNE. 

5 Summary 

A calculation of the multiplicity of heavy quarks from gluon splitting in e^e~ annihilation, as 
a function of the heavy jet mass, has been presented. The shape of the result is similar to that 
predicted by Monte Carlo event generators at the Z*^, lying between the different models, but in 
better agreement with ARIADNE at higher energy. 
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